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Abstract 

A model of competition between complete fusion and quasifission channels 
in fusion of two massive nuclei is extended to include the influence of dissi- 
pative effects on the dynamics of nuclear fusion. By using the multidimen- 
sional Kramers-type stationary solution of the Fokker-Planck equation, the 
fusion rate through the inner fusion barrier in mass asymmetry is studied. 
Fusion probabilities in symmetric 90 Zr+ 90 Zr, 100 Mo+ 100 Mo, 110 Pd+ 110 Pd, 
136 Xe+ 136 Xe, almost symmetric 86 Kr+ 136 Xe and 110 Pd+ 136 Xe reactions are 
calculated. An estimation of the fusion probabilities is given for asymmetrical 

62 M+ 208p b) 70 Zn+ 208p b) 82g e+ 208p b) and 48 Ca+ 244p u reactions used for the 

synthesis of new superheavy elements. 
PACS:25.70.Jj, 24.10.-i, 24.60.-k 
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I. INTRODUCTION 



The competition between the complete fusion and quasifission processes occurs in the 
reactions with massive nuclei at bombarding energies smaller than 15 MeV/nucleon. In these 
reactions the quasifission channel dominates and leads to a strong reduction of few orders 
of magnitude of the fusion cross section a C N [0]- The competition between the complete 
fusion and quasifission processes was not considered in the calculations with various models, 
namely with the macroscopic dynamical 0, optical || and surface friction [|J models. It 
was shown in |TJ that these models do not lead to correct fusion cross sections for reactions 
with heavy nuclei The new model suggested in Jy] yields a good agreement between 

the theoretical predictions and experimental data on the fusion of heavy nuclei. Within this 
model the fusion cross section acN is defined as 

Jmax 

&CN{Ec. m ) = X! a c{E c .m., J)PcN{E c . m ., J). (1) 
J=0 

Here, a c is the capture cross section at the bombarding energy E c . m . and angular quantum 
number J and Pcn the fusion probability after the capture stage of collision, which takes the 
competition between the fusion and quasifission processes into account. The value of J ma x 
depends on E cia _ and is smaller than Jb } =o & t which the fission barrier in the compound 
nucleus vanishes ||. In order to calculate a c , either optical || or surface friction || models 
can be used Q. Pcn = 1 is supposed in the models of Refs. HH]. In the macroscopic 
dynamical model Pcn = 1 for E c _ m , > Be + E xx and Pcn = for E c , m _ < Be + E^ 
where Be and E xx are the Coulomb barrier and the extra-extra push energy, respectively. 
The main advantage of the model fl[] in contrast to the models is the treatment of 

the competition between the complete fusion and quasifission in the calculation of Pcn- In 
the present paper the model |]] is extended by including dissipative effects on the fusion 
dynamics. 

The fusion process is considered in |JJ as the evolution of the dinuclear system (DNS) in 
which nucleons are transfered from the light nucleus to the heavy one . The DNS-concept 
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is based on the information obtained from the investigation of deep inelastic collisions of 
heavy ions |TU| . The DNS is formed at the initial stage of the reaction when the kinetic 
energy is transformed into the excitation energy of the nuclei. The initial DNS is localized 
in the minimum of the pocket of the nucleus- nucleus potential V(R) (Fig. [I]) at R — R m 
where R is the relative distance between the interacting nuclei. In our approach we have 
R m ~ -R1 + -R2 + O.5 fm where R\ and R2 are the radii of the nuclei in the DNS. Then the DNS 
evolves by a diffusion process in the mass asymmetry degree of freedom 77 = (A\ — A<i)jA 
to the compound nucleus and fuses (Fig. [I]). A\ and A2 are the mass numbers of the nuclei 
and A = Ai + A2. Besides the motion in 77 a diffusion process in the relative distance occurs. 
This process leads to the decay of the DNS which we denote as quasifission. For quasifission 
the DNS should overcome the potential barrier {B q f) which coincides with the depth of the 
pocket in the potential V(R) (Fig. [I]). The important peculiarity of the DNS evolution to the 
compound nucleus is the appearance of an inner fusion barrier Bj us in the mass asymmetry 
degree of freedom with its top (the Businaro-Gallone point) at 77 = tjbg which coincides with 
the maximum of the DNS potential energy as a function of 77 (Fig. [TJ). The value of Bj us 
supplies the hindrance for the complete fusion in the DNS-concept. The energy required to 
overcome the fusion barrier Bj us is contained in the DNS excitation energy. As was found in 
|l|.[Tl|Jl^] , the value of Bj us can be much smaller than the extra-extra push energy predicted 



m 0- 

The dissipative large-amplitude collective nuclear motions, which occur in fission, quasi- 
fission, fusion and heavy-ion reactions, can be analyzed within the transport theory. It is 
known that the Kramers expression |14J] yields an excellent approximation to the fission 



rate in comparison to the solution of the Fokker-Planck and Langevin equations |[T5| , pO| , p!6 
The fusion probability Pqn (|3) can be obtained with the Kramers-type expression for the 
rate of fusion through the inner fusion barrier. The main advantage of the Kramers-type 
expression for the fusion rate in contrast to the statistical method is the possibility to 
include nuclear viscosity in fusion process. In this paper, the two-dimensional (R and 77) 
Kramers- type expression is applied to the calculation of Pcn- In the following we analyt- 



ically demonstrate the dependence of the fusion probability on the values of the friction 
coefficients, the temperature T of the DNS, and the values of Bj us and B q f. As examples, 
the symmetrical, almost symmetrical, and asymmetrical reactions will be considered. 



II. MODEL 



A. Fusion probability 



The variables R and i] are the relevant collective variables used in [|I],|T7j] to describe the 



DNS evolution. The neck degree of freedom P, |l8| , |l9f which is important in macroscopic 



dynamical models is not considered in this paper. As follows from our analysis |1J, the 



neck is not a relevant collective variable for values of R > Ri + R2 which are important in 
the DNS. The size of the neck in the liquid drop model and the overlap region in the frozen 
density approximation are close to each other for R\ + R2 < R < -Ri + R2 + 1 fm [|19[ . Based 
on this fact and the time scale, which is of interest, we assume that the individuality of the 
DNS nuclei in the fusion process is retained, and follow the approach suggested in |P],|I2| . 

Since the fusion probability increases with the mass asymmetry, the DNS with larger 
mass asymmetry supplies favorable conditions for complete fusion. We assume that the 
fusion occurs inevitably for 77 > i]bg- The diffusion in 77 is important in our consideration 
of the fusion process. The simultaneous investigation of the transport processes in 77 and R 



variables allows us to calculate the fusion and quasifission probabilities | 12| . In our opinion, 
the fusion probability is small in many reactions with heavy nuclei because the initial mass 
asymmetry r\i is smaller than tjbg and the diffusion in rj is much smaller than the diffusion 
in R. 

The leakage of probability through the fusion barrier in rj is defined by the rate X v (t) at 
V — Vbg (Fig- Then we obtain 

Pcn = i \(t)dt (2) 







Here, t n is the life-time of the DNS with 



to 

J[\ R (t) + \ v (t)]dt = l. (3) 



For symmetrical and almost symmetrical systems, the symmetry of the fusion process with 
respect to rj = should be taken into account. Then the initial DNS is near the minimum 
of the potential energy U{rj) which is a symmetric function with respect to rj = (Fig. p. 
Fusion occurs when DNS reaches the barriers at either rj = t]bg or rj = — i]bg- The time 
dependence of the rates A*(t) (i = R,rj) can be taken in the following way 

Ht) = Af r (^zt^(^ + o(t - tj)) , (4) 

where Xf r are asymptotic values of the fusion or quasifission rate Xi(t) at the corresponding 
barriers (Fig. [l|) and 9{t) is a step function. Here, we assume that after the exponential 
growth during the transient time Tj the rate \{t) reaches the asymptotic value. 
Using (f|) we obtain from Eqs. (0) and (|3]) 

\Kr \Kr\Kr 
p _ pO A p _ \ \ A R Ty ~ Tr 



1 \% r T R + \« r T v 

where (3 = e — 1 « 1.72. Assuming instead of (^) the linear grows of A,(t) |]TB| 

A i (t)=Af r (^(r < -t)+e(t-r i )) (7) 

we obtain the Eqs. (j5|) and (|6|) but with (3 = 2. The first terms in (|5]) and (||) are the 
contributions from the quasistationary rates. The second terms are related to the transient 
time. One can see from fl5|) that we can neglect the transient time for r» <C 1/Af r (i = i?, 77) 
or tr ~ r^. This is true for all reactions under consideration excepting the 136 Xe+ 136 Xe 
and 110 Pd+ 136 Xe reactions where the differences between Bj us and B q f are very large. Note 
that the role of the transient stage decreases with decreasing excitation energy of the DNS 
because the exponential increase of 1/Afj r is larger than the logarithmic increase of the 
transient times (see equations below). 



To obtain the asymptotic fusion and quasifission rates A 



Kr 



R,T]), we use the for- 



malism elaborated in Refs. p~5|JT6[1 . We approximate the expression for the quasistationary 
rate \f r over a multidimensional potential barrier |16| with a Kramers-type formula 



A 



Kr 



(r/ft) 2 

, .Br B v 



+ 4- 



(r/ft) 2 

, .Br B v 



v. 








exp 


T. 



(8) 



Here, Bi (i = R,rj) defines the height of the fusion (B v = B*j us ) or quasifission (B R = 
B q f) barriers. The possibility to apply the Kramers-type expression to relatively small 



barriers (Bi/T > 0.5) was demonstrated in [20|. The local thermodynamic temperature 
T is calculated with the expression T = 



<E* /a, where a = A/12 MeV -1 and E* is the 
DNS excitation energy. In Eq. (]8p, the frequencies cu i J (j = R,r)) of the inverted harmonic 
oscillators approximate the potential in the variables i = R, 77 on the tops of the barriers Bj, 
and Ui are the frequencies of the harmonic oscillators approximating the potential in the same 
variables for the initial DNS. Since the local oscillator approximation of the potential energy 
surface is good for the reaction considered, we neglected the nondiagonal components of the 
curvature tensors in (||). In our calculations, we use the simple approximate expressions for 
the friction coefficients 

7«/ = Yn w /K, = R,rj), 

which were obtained by the linear response theory pi ]. The quantity V denotes an average 
double width of the single-particle states. The calculation of the mass parameters hrr and 
fi vv is given in [|T2| , p2| . The role of the nondiagonal components of the tensors of inertia 



and friction depends very much on the choice of the collective variables | 21]j22 |. By using 
the collective variables R and 77 for describing the evolution of the DNS we neglect the 
nondiagonal component of the tensor of inertia because hr v <^ ^/JlRRprrn f° r \v\ < \Vbg\ 
. The nondiagonal mass coefficient fj,R V = in the DNS takes an essential role only 



for \7)\ > \t]bg\- As was shown in |T^,|T9[, the friction coefficients 7^ and 7 W obtained 
with T = 2 MeV have the same order of magnitude as the ones calculated within the other 
approaches. 
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B. Potential energy of the DNS 



The value of 



0J4 



\ 



d 2 U(R,r], J) 



di 2 

is easily calculated with the potential energy of the DNS 



U(R, V , J) = B 1 + B 2 + V(R, J) - [B 12 + V; ot (J)}- (9) 

Here, B\, B 2 , and B\ 2 are the binding energies of the fragments and the compound nucleus 
and are calculated with liquid-drop masses for large excitation energies and with realistic 
masses |23] for small excitation energies. The isotopic composition of the nuclei forming 
the DNS is chosen with the condition of a iV/Z-equili brium in the system. The value of 
U (R, r], J) is normalized to the energy of the rotating compound nucleus by Bi 2 + V' rot . The 
nucleus- nucleus potential V(R, J) in (^) is calculated as described in |23J. The retaining 
individuality of the DNS nuclei during the time which is of interest allows us to calculate 
the DNS potential energy by the method presented in |24j] . The calculated driving potential 



U(R m ,r], J = 0) = U(rj) as a function of r\ and the nucleus-nucleus potential V(R, J = 0) 
as a function of R in the reaction 90 Zr + 90 Zr is presented in Fig. |l[ 

The potential energy of the DNS as a function of 77 and R depends on the temperature, 
shell effects and angular momentum. In the present paper we do not analyse in details 
the dependence of the potential energy surface on temperature because this demands the 
introduction of an additional parameter. We distinguish two cases: The first case corre- 
sponds to large excitation energies of the initial DNS when the liquid drop binding energies 
and spherical shapes of the nuclei in the DNS can be used in the calculations. The second 
case corresponds to the cold fusion with small E* when the realistic binding energies are 
taken in (H). For small values of E*, the deformation of the nuclei in their ground states 



25fl is taken in the DNS into account when the barrier heights are calculated. In order to 
demonstrate the influence of the shell and deformation effects at small E*, the calculation 
of U(R m ,r], J = 0) is presented in Fig. ^| for the 86 Kr+ 136 Xe reaction. In this reaction the 



deformation of the nuclei in the DNS leads to a decrease of Bj us . After the smoothing over 
small oscillations caused by even-odd effects we can use the expression @ for Af r in first 
approximation. The calculated driving potentials for the reactions leading to superheavy 
nuclei are presented in Below we discuss the case of the 86 Kr+ 136 Xe reaction where 
the fusion probability is calculated with driving potentials based on spherical and deformed 
shapes of the nuclei. The calculations of complete fusion probabilities with spherical (the 
liquid drop binding energies) and deformed (the realistic binding energies) nuclei give us 
the interval where realistic values of Pcn are situated. As noted in Ref. [|l^], the values of 
Pcn for small E* can be smaller or larger than the ones for large E*. In our calculations 
of the fusion probability, the following values hu^ R ~ 0.8 — 1 MeV, hu^ 1 ~ 3 — 3.5 MeV, 
hwjj* w 1 - 1.5 MeV, hu^ « 1.5 - 2 MeV, hu R w 1.5 - 2 MeV, and hu v « 0.8 - 1 MeV 
are used for the reactions considered. 

The dependences of the barrier heights on J for the reactions 90 Zr+ 90 Zr and 110 Pd+ 110 Pd 
are presented in Fig. []. The values of Bj us and B q f are slightly changed when J in- 
creases from to 25ft. For heavier systems, this changes are certainly smaller because of 
the larger moment of inertia. As a result, the values of Pcn slightly differ from those 
calculated with J = 0. In order to calculate a c (E cm ,J) in ([!]) the simple expression 
&c{E c .m.i J) — 7rA 2 (2J+ l)T(i? c . m ., J) is often used (A is the reduced de Broglie wavelength). 
The transmission coefficient T(E c m , J) through the Coulomb barrier restricts also the an- 
gular momentum range. The weight function (2J + 1)P C n{E cu11 J) used to obtain ctcn 
has its maximum at J = 20h. Since the calculation of (Jcn is of interest to determine the 
evaporation residues cross sections, only low angular momenta can be considered. Indeed, 
the surviving probabilities of the compound nuclei in the reactions considered are narrow 
functions peaking at all energies at J values in the vicinity of zero |5| . Although the precise 
calculation of the excitation function demands the dependence of Pcn on J, for the esti- 
mations of the evaporation residues cross sections, we can use the values of ctcn calculated 
with Jmax = 10 — 15fr and PcN{E c . m ., J = 0) with a good accuracy. 



S 



III. RESULTS AND DISCUSSIONS 



A. Effect of transient times 

Since the difference between Bj us and B q f and their absolute values can be large in 
some reactions, the role of the transient time in the calculation of Pcn and to should be 
estimated. The motions in R and rj are close to the underdamped and overdamped motions, 
respectively. Therefore, the transient time for the realistic values of T can be estimated by 



the expressions fl3 



T, 



r R = ^\n(lOB qf /T), (10) 



ln(10B; M /T), (11) 



The calculated time dependences \ n {t) with Eq. (24) in |L2| are presented in Fig. f| for the 
reactions 90 Zr+ 90 Zr and 136 Xe+ 136 Xe at Y = 2 MeV and E* = 30 MeV. We find that the 
transient times taken from Fig. [| are practically the same as the ones calculated with (|ITD. 
The calculated transient times and life-times for different systems are given in Table 1. One 
can see that the transient stage does practically not effect the values of to- 

It is known [27| that the consideration of the transient stage is important for large 



excitation energies when the particle emission change the system during the fission or fusion 
time. However, this is not so in our cases and the calculation with the transient time 
(Table 2) leads only to a decrease of Pcn by maximal 30% in the reactions 136 Xe+ 136 Xe 
and 110 Pd+ 136 Xe in comparison with the value of P CN in (||) for T = 2 MeV. In other 
reactions considered the effect of the transient stage is negligible. Therefore, we can neglect 
the transient time in the calculations of Pcn for the most considered reactions. For the 
reactions with Bj us —B q f > 15 MeV, the values of Pcn are small and the corrections APcn 
from the transient stage are practically within the accuracy of the calculated barrier heights. 
In these reactions the values of Pcn calculated without APcn in (|5|) are the upper limits of 
the fusion probabilities and can also be used for estimations of fusion cross sections. 
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B. Large excitation energies of the DNS 



The values of Pcn (Table 3) are in agreement for most of the reactions with the ones 
extracted from the experimental data 17. 2^.29] and with the results of our previous calcula- 
tions HUTTJ. For example, for the reactions 90 Zr+ 90 Zr, 100 Mo+ 100 Mo, and 110 Pd+ 110 Pd the 
values of Pcn ~ 4 x 10 _1 , 10~ 2 , 10~ 4 , respectively, yield a good agreement with the experi- 
mental data of crcN{E c .m.)- For the 110 Pd+ 110 Pd reaction, the value of <JcN{E c . m .) calculated 
with the macroscopic dynamical model M is about three orders of magnitude larger than 
the experimental one. Thus, the competition between the complete fusion and quasifission 
processes is extremely important in the DNS evolution. Perhaps, due to the small value of 
Pcn in the 110 Pd+ 136 Xe reaction, the fusion was not observed in [p9fl . 

The values of Pcn can be also calculated from @ by using the one-dimensional Kramers- 
type expression instead of Eq. (§) : 

1 <jji 



A 



Kr 



2tt uf 1 



^) +(^) 2 -^) ex P 



Bi 
' T 



(12) 

Therefore, the 



We approximately find the same results with this formula as with Eq. (|_ 
estimations of the transient times with fliTf ) and ([11]) are realistic. 

In addition to the reactions presented in Table 3 the fusion probabilities were calculated 
for the reactions 96 Zr+ 124 Sn and 124 Sn+ 124 Sn with T = 2 MeV, J = and the DNS excitation 
energy 30 MeV. The values of Pcn as a function of Z\ x Z 2 {Z\ and Z2 are the charge numbers 
of the colliding nuclei) are presented in Fig. ||]. One can see the exponential decrease of the 
fusion probability with increasing Z\ x Z2 in the symmetric and almost symmetric reactions. 
Therefore, the experimentally observed || rapid fall-off of the fusion cross sections with 
increasing Z\ x Z 2 is simply explained in our model. 

The fusion and quasifission rates decrease with increasing T. However, Pcn increases 
quickly (the quasifission rate decreases more strongly than the fusion rate) and reaches a 
plateau at T m 4MeV because the changes of the dissipative effects in 77 and R variables start 
to compensate each other. The calculated dependence of X^ r , Afj r and Pcn on the friction 
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parameter T is shown in Fig. |^ for the 110 Pd+ 110 Pd reaction. Since the fusion probability 
increases not much with T, our calculations agree with the values of Pcn obtained within 
the approach [[L]] based on statistical assumptions within an order of magnitude. It is seen 
that the results of the calculations are not crucial to the exact value of the parameter Y for 
T > 2 MeV. A realistic value of the parameter Y is about 2 MeV. This value should be used 
in the further calculations. 

The dependences of Pcn on the excitation energy E* = E c , m , — V(R m ) of the DNS 
are given for the 110 Pd+ 110 Pd and 86 Kr+ 136 Xe reactions in Fig. 0. We observe that Pcn 
increases with E* because the value of X^ r in (Q) increases slower than X^ r grows. For the 
86 Kr+ 136 Xe reaction, the Pcn values were calculated for two cases. In the first case, liquid 
drop masses in (§) and spherical nuclei in the DNS were used to calculate Bj us and B q f. In 
the second case, realistic masses in (|9|), spherical nuclei in the initial DNS and a deformed 
heavy nucleus in the ground state near r\ = t)bg were taken. The pole orientation of the nuclei 
leads to the minimum of the potential energy in the DNS. Due to the deformation effect near 
V — Vbgi the value of V(R m ) decreases as compared to the calculation with the spherical 
nuclei and U(R m ,rj BG ) decreases. As a result, BJ US decreases and Pcn increases (Fig. ^). 
The driving potential with realistic binding energies and deformation effects is preferable for 
small excitation energies. The driving potential with the liquid drop binding energies is good 
for large excitation energies. For the 5n channel in the 86 Kr+ 136 Xe reaction, the excitation 
energy of the compound nucleus is about 46 MeV (E* = 30 MeV) and P C n = 4 x 10~ 2 . With 
a c = 23 mb estimated with the model || and average value of < Y n /Yf >= 0.3 taken from 
Ref. |3(J, we obtain the evaporation residue cross section o ER « <j c Pcn < r„/Ty > 5 = 2.2 /ib 



which is in agreement with the experimental value 5 fib [^T 



The energy threshold for complete fusion is related to the fusion barrier B*j us (see Table 3) 
and can be much smaller than the extra-extra push energy which, for example, is E^= 60 
MeV and 30 MeV in 110 Pd+ 110 Pd and 62 Ni+ 208 Pb reactions, respectively, predicted in the 
macroscopic dynamical model 0. This result of our model is in agreement with recent 
experimental data on the synthesis of the new superheavy elements |5|,§]. 
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C. Small excitation energies of the DNS 



For low excitation energies of the initial DNS in asymmetrical reactions 62 Ni+ 208 Pb 
^ 270 110, 70 Zn+ 208 Pb ^ 278 112, 82 Se+ 208 Pb ^ 290 116, and 48 Ca+ 244 Pu ^ 292 114 which 
were used for the synthesis of new superheavy elements, the initial DNS is in the local 
minimum of U(rj) due to shell effects (realistic binding energies are taken in (01)) ||11|| . In this 
case, we can use the Kramers-type expressions (P|) to estimate the values of Pcn (Table 4). 
For these reactions, deformations of the DNS nuclei were taken into account • Since 

in the 48 Ca+ 244 Pu reaction the heavy nucleus is deformed even in the initial DNS, the 
treatment of the deformation of nuclei leads to larger values of Bj us as compared to the 
calculation with spherical nuclei (Table 4). Deformation effects lead to a decrease of Bj us in 
other reactions. The detailed discussion of the influence of the deformation and orientation 



of the DNS nuclei on the value of Bj us is given in f!2 |. In the reactions leading to the 
superheavy nuclei only partial waves with very small J up to J max = 10 — 15h contribute 
the evaporation residue cross section because these nuclei are instable against fission for 
larger J. The effect of the transient time in these reactions seems to be very small. 

Using the results presented in Table 4 one may explain the smaller fusion yields of the 
nuclei with Z — 112 as compared to the fusion yields of the nuclei with Z = 110 ||. As given 
in Table 4, the probability to obtain a superheavy nucleus with Z — 116 in the 82 Se+ 208 Pb 
reaction is very small. The use of this combination for producing heavy compound nuclei 
with a small excitation energy may be problematic. In spite of the larger value of Pcn in 
the 48 Ca+ 244 Pu reaction as compared with others in Table 4, the compound nucleus seems 
to be more excited due to the Q-value. The analysis of the surviving probability W sur of the 
compound nucleus is extremely important to estimate the yield of the element with Z — 114 
in this reaction. 

By comparing Pcn in Tables 3 and 4 with Pcn given in the experimental work, we 
should bear in mind that the experimental values were extracted from the evaporation 
residues cross sections cr ER (E cm ) by model assumptions about the surviving probability 
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W sur of the excited compound nucleus (iter = &cPcNW sur ) |I|,p8|. For the In channel in 
the 62 Ni+ 208 Pb reaction, the excitation energy of the compound nucleus is about 13 MeV 
and Pcn = 7 x 10~ 6 . Extrapolating the systematic representation in Refs. p0|,|32| for the 



nucleus 270 110 we estimated W SUT ~ T n /Tf = 3 x 10 -4 . With a c = 4 mb estimated with 
the optical model |3j and the values of Pcn and W mr , we obtain o ER = 8.4 pb which is in 
agreement with the experiment ||. 

A discussion on the accuracy of the model presented is necessary for applying it to 
reactions producing superheavy elements with very small cross sections. As in any model, 
certain assumptions are used in our approach. However, with the same assumptions and 
the set of the parameters our model is able to describe the experimental data for different 
reactions. As it is seen, our model gives good results for the 90 Zr+ 90 Zr reaction for which the 
known traditional models are applicable. However, our model describes also the experimental 
data well in the case of reactions where the fusion cross sections are very small and other 
models may fail [|TJ. 

IV. SUMMARY 

The new model suggested to calculate the probability of the fusion of heavy nuclei is 
useful for the analysis of experimental data. The results obtained support the use of the 
simple statistical assumptions of our previous studies p]Jll|. A good description of Pcn 



in our model can be considered as an evidence for the DNS-concept providing a realistic 
interpretation of the mechanism of fusion process. In our opinion the fusion probability is 
small in many reactions with heavy nuclei because quasifission plays a major role. Without 
taking the quasifission into account, the explanation of the experiments on the fusion of heavy 
nuclei is not possible. Based on the results presented, we plan calculations of the evaporation 
residues cross sections for a large set of reactions used to produce the superheavy elements. 
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TABLES 

TABLE I. Calculated transient times (|lO|) and (11), quasistationary values of 1/Af r and life 
time (|6|) in the symmetric and almost symmetric reactions for J = 0, L = 2 MeV and the excitation 
energy E* = 30 MeV of the initial DNS. The calculations were made with Bj us and B q f given in 
Table 3. 



Reactions 


TR 


i/xf- 

' it 


'1 


< V 


^00 


At 




xl0~ 21 s 


xl0~ 21 s 


xl0- 21 s 


xl0" 19 s 


xl0~ 21 s 


xl0~ 21 s 


90 Zr+ 90 Zr 


1.2 


104 


2.4 


1.6 


63 


1.0 


ioo Mo+ ioo Mo 


1.1 


60 


2.8 


39 


58 


0.6 


iio pd+ iio pd 


1.0 


32 


3.1 


2800 


32 


0.6 


86 Kr+ 136 Xe 


1.1 


70 


2.8 


18 


68 


0.6 


110p d+ 136 Xe 


0.5 


50 


3.2 


8300 


5 


0.3 


136 Xe+ 136 Xe 


0.5 


50 


3.5 


7.1xl0 6 


5 


0.3 



TABLE II. Calculated values of P$ N , AP CN and the ratio AP cn /Pcn ( see Ec l- ©) in the 



symmetric and almost symmetric reactions for J = 0, L = 2 MeV and the excitation energy 
E* = 30 MeV of the initial DNS. The calculations were made with B% us and B q f given in Table 3. 



Reactions 


pO 

CN 


AP C N 


AP C n/P Q cn 


90 Zr+ 90 Zr 


4.0 x 10" 1 


2.7 x 10" 3 


6.7 x 10~ 3 


ioo Mo+ ioo Mo 


1.5 x 10~ 2 


2.6 x 10~ 4 


1.7 x 10~ 2 


iio pd+ iio pd 


1.1 x 10" 4 


4.2 x 10~ 6 


3.8 x 10~ 2 


86 Kr+ 136 Xe 


4.0 x 10~ 2 


5.6 x 10" 4 


1.4 x 10~ 2 


110p d+ 136 Xe 


5.5 x 10~ 6 


1.7 x 10" 6 


3.1 x 10" 1 


136 Xe+ 136 Xe 


6.1 x 10~ 9 


1.9 x 10" 9 


3.1 x 10" 1 
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TABLE III. Calculated fusion probability Pcn i n t ne symmetric and almost symmetric reac- 
tions for different friction parameters V. The calculations were made for J = and the excitation 
energy E* = 30 MeV of the initial DNS. For the calculated values of Pcn , liquid drop masses and 
spherical nuclei in the DNS were used in @. 



Reactions 


R* 
B fus 

^IVlt; V J 


\ ivie v i 


Pcn 


T=0 MeV 


T=l MeV 


T=2 MeV 


T=3 MeV 


T=4 MeV 


90 Zr+ 90 Zr 


6 


5 


3-lxlO" 1 


3.6X10- 1 


4.0X10" 1 


4.5X10" 1 


4.8xl0 _1 


ioo Mo+ ioo Mo 


10 


4 


8.3xl0~ 3 


l.OxlO" 2 


1.5xl0- 2 


1.8xl0~ 2 


2.0xl0~ 2 


iio pd+ iio Pd 


15 


3 


6.5xl0~ 5 


7.9xl0" 5 


l.lxlO" 4 


1.4xl0~ 4 


1.5xl0~ 4 


86 Kr+ 136 Xe 


8.5 


4 


2.3xl0~ 2 


2.7xl0~ 2 


4.0xl0~ 2 


5.0xl0~ 2 


5.6xl0~ 2 


110p d+ 136 Xe 


15.5 


0.5 


9.0xl0~ 7 


2.2xl0~ 6 


3.8xl0~ 6 


4.4xl0~ 6 


5.0xl0~ 6 


136 Xe+ 136 Xe 


22.5 


0.5 


1.6xl0~ 9 


2.5xl0~ 9 


4.2xl0~ 9 


4.9xl0~ 9 


5.7xl0~ 9 
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TABLE IV. The same as in Table 3, but for asymmetric reactions used for producing super- 
heavy nuclei. The deformation effects are taken into account (see text). The calculations were 
made for J = and the excitation energy E* = 15 MeV of the initial DNS with exception of the 



Ca+ Pu reaction. For the Ca+ Pu reaction, the choice of the initial excitation energies 
in the calculations of Pcn with (E* = 33 MeV) and without (E* = 15 MeV) deformation effects 
(sph.) yields the same excitation energy of the compound nucleus (« 40 MeV). 



Reactions 


&fus 

(MeV) 


B qf 

(MeV) 


Pcn 


T=0 MeV 


T=l MeV 


T=2 MeV 


T=3 MeV 


T=4 MeV 


62 Ni+ 208p b _^270 UQ 


8 


1.5 


1.4xl0~ 4 


1.6xl0~ 4 


2.2xl0~ 4 


2.9xl0~ 4 


3.2xl0~ 4 


70 Zn+ 208p b ^278 m 


9.5 


1 


9.1xl0~ 6 


l.lxlO- 5 


1.5xl0- 5 


1.9xl0- 5 


2.2xl0~ 5 


82 Se+ 208p b ^290 U g 


12.5 


0.5 


9.4xl0~ 8 
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2.2xl0~ 7 


48 Ca+ 244p u ^292 m 


12 
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8.0xl0~ 4 


9.1xl0~ 4 


(sph.) 48 Ca+ 244 Pu ^ 292 114 
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2.3xl0~ 3 


2.8xl0~ 3 


3.9xl0~ 3 


5.0xl0~ 3 


5.6xl0~ 3 
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FIGURES 

FIG. 1. Calculated potential energy of the DNS corresponding to the 180 Hg compound nucleus 
as a function of the mass asymmetry for zero angular momentum (upper part). The energy scales 
are normalized to the total energy of the compound nucleus. Liquid drop binding energies are used 
in (^). Calculated nucleus-nucleus potential in the 90 Zr+ 90 Zr reaction for zero angular momentum 
(bottom part). 

FIG. 2. Calculated potential energy of the DNS in the 86 Kr+ 136 Xe reactions as a function of 
r\ for J = 0. The result calculated with liquid drop binding energies in @ is presented by the 
dashed line. The driving potentials with and without deformation effects in the DNS, obtained with 
realistic binding energies in @, are presented by thick and thin lines, respectively. The deformed 
nuclei of the DNS are assumed in the pole orientation. 

FIG. 3. Dependences of B qf , B* fus , P CN and (2J + l)P C N on J for the reactions 90 Zr+ 90 Zr 
(solid lines) and 110 Pd+ 110 Pd (dashed lines). In the calculations of Pcn the values E* = 30 MeV 
and r = 2 MeV are used. 

FIG. 4. Calculated fusion rates as a function of time for the reactions 90 Zr+ 90 Zr (solid line) 
and 136 Xe+ 136 Xe (dashed line) at E* = 30 MeV and T = 2 MeV. 

FIG. 5. Calculated values (solid points) of the fusion probability as a function of Z\ x Z<i for 
the reactions presented in Table 3 and the reactions 96 Zr+ 124 Sn (Bj us = 9 MeV, B q j = 4 MeV) 
and 124 Sn+ 124 Sn (B* fus = 16 MeV, B qf = 0.5 MeV) at J = 0, T = 2 MeV and E* = 30 MeV. Z x 
and Z2 are the charge numbers of the colliding nuclei. The solid line is drawn to guide the eye. 

FIG. 6. Fusion (long-dashed line) and quasifission (short-dashed line) rates Af r (i = r],R) and 
fusion probability Pcn (solid line) as a function of the friction parameter T for the 110 Pd+ 110 Pd 
reaction at J = 0. 
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FIG. 7. Fusion probability Pcn as a function of the excitation energy E* = E c , m , — V(R m ) 
of the initial DNS for the reactions 110 Pd+ 110 Pd (upper part) and 86 Kr+ 136 Xe (bottom part) at 
J = 0. The results were obtained by using liquid-drop masses and spherical nuclei (solid lines) and 
realistic masses and a deformed heavy nucleus (dashed line) in the DNS. The calculations were 
made with a friction parameter T = 2 MeV. 
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Fig. 2 NUCPHA 1629 
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Fig. 3 NUCPHA 1629 
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Fig.4 NUCPHA 1629 
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Fig. 5. NUCPHA 1629 



27 




2 4 6 

T (MeV) 



28 




10 20* 30 40 50 

E (MeV) 



29 



